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A spatially flat Robertson- Walker spacetime driven by a cosmoiogicai constant is non-conformally 
coupled to a massless scalar field. The equations of semiclassical gravity are explicitly solved for this 
, case, and a self-consistent de Sitter solution associated with the Bunch-Davies vacuum state is found 

. (the effect of the quantum field is to shift slightly the effective cosmoiogicai constant). Furthermore, 

' it is shown that the corrected de Sitter spacetime is stable under spatially-isotropic perturbations 

of the metric and the quantum state. These results are independent of the free renormalization 
parameters. 
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I. INTRODUCTION 



I Our present understandiiig of cosmology assumes that the universe underwent a short period of accelerated ex- 

O , pansion known as inflation [il, i, i, i, H • The inflationary scenario has been remarkably successful in explaining the 
observed anisotropics of the cosmic microwave background [1, 0, H [1] ■ In most inflationary models the accelerated 
^ ] expansion phase is close to but never exactly de Sitter and this phase eventually ends when the kinetic energy of 
the inflaton field driving inflation starts to dominate over the potential term. On the other hand, observations of 
distant supernovae indicate that the universe is presently under goin g a period of accelerated expansion [Tol . [Tl| that 
,—1 , may be driven by a small non- vanishing cosmoiogicai constant l^, H, [3j If that is the case, the geometry of 



J> ■ our universe would tend to that of de Sitter spacetime at sufficiently late times. Thus, a detailed knowledge of the 
CN [ physics associated with de Sitter space may play a key role in understanding both the very early universe as well as its 
ultimate fate. Furthermore, it is conceivable that studying a possible screening of the cosmoiogicai constant driving 
de Sitter space, due to quantum effects, could shed some light on the huge fine-tuning problem that the current value 
. of the cosmoiogicai constant seems to pose. 

An open question which has recently received increasing attention is whether the quantum fluctuations of the metric 
and the matter fields in de Sitter space can give rise to large back-reaction effects on the mean background geometry. 
It has been argued that in pure gravity with a cosmoiogicai constant the infrared effects due to two graviton loops 
and higher-order radiative corrections could lead to a secular screening of the cosmoiogicai constant [H, |T3| ■ There 
have also been proposals that a significant screening of the cosmoiogicai constant could a ppe ar in chaotic inflationary 
models at one loop when both the metric and inflaton field fluctuations are considered [iSl [l9l . [20l . [2l| . In all these 
cases the quantum fluctuations of the metric play an essential role. However, whenever the metric perturbations are 
' quantized, one needs to confront the problem of defining proper diffeomorphism-invariant observables in quantum 
gravity [22l |. even when treated as a low-energy effective field theory, fn particular one needs to make sure that 
the secular screening found in the analysis mentioned above is not simply a gauge artifact. As a matter of fact, 
it was shown in Refs. H^l that when a suitable gauge- invariant measure of the expansion rate was considered 
the screening effect previously found in chaotic infiationary models was not actually present (at least for single field 
models). Similarly, a recent reanalysis of the pure gravity case which made use of a diffeomorphism-invariant measure 
of the change of the expansion rate revealed the absence of secular effects to all orders in perturbation theory ^ . 

In recent work it was also found that the back reaction due to one-loop effects of massless non-conformal fields 
can give rise to substantial deviations from de Sitter spacetime [1^ [13] . This would be the case even if the quantum 
fiuctuations of the metric are not considered (provided that there is some other massless non-conformal field in 
addition to gravitons) . Not quantizing the metric perturbations means that the ambiguity associated with the gauge- 
fixing term for the metric perturbations is no longer present, and the mean geometry is a perfectly well-defined (and 
gauge-independent) object. Hence, the subtleties mentioned in the previous paragraph do not apply. The heuristic 
argument provided in Refs. [2^ . [27| is that one-loop contributions from massless non-conformal fields correspond to 
logarithmic non-local terms (conformal fields only produce local terms for Robertson- Walker geometries), and by 
analogy with the situation in pion physics, one expects that they become important in the infrared limit. However, 
it is not clear that pion physics constitutes a good analogy because the derivative coupling of the matter fields to the 
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metric generates higher powers of the momentum. This point can be illustrated with the simple example of small 
metric perturbations around flat space. In that case the Fourier transform of the inverse propagator behaves like 
+ 6 (/c^/mp) ln(fc^//iQ)], where & is a dimensionless number roughly of order one, nip is the Planck mass and /zq 
is some fixed mass scale. One can see that although the logarithm grows in the infrared, the whole term actually 
decreases because it is suppressed by the factor (fc^/rrip). Indeed, a detailed calculation of the quantum radiative 
corrections to the Newtonian potential shows that the contribution from that term is suppressed by the square of the 
ratio of the Planck length over the radial distance [2^[2^,[33]- Of course the case of a Robertson- Walker (RW) metric, 
which involves a time-dependent scale factor, is not so simple and deserves a careful analysis in order to compare with 
the detailed calculation in Refs. [26l. [27j. 

In this paper we approach this problem by explicitly solving the back reaction on the mean gravitational field 
due to the quantum effects of a massless non-conformally coupled scalar field when the quantum fluctuations of the 
metric are not considered. This kind of one-loop calculation is entirely equivalent to studying the corresponding 
back-reaction problem in the semiclassical gravity framework jsT], [H, [s^ by solving self-consistently the semiclassical 
Einstein equation, which includes the suitably renormalized quantum expectation value of the stress tensor operator 
acting as a source. Specifically, in our calculation we assume the presence of a cosmological constant, which would 
lead to a de Sitter solution in the absence of quantum effects, and simplify the problem by focusing on RW geometries, 
corresponding to spatially homogenous and isotropic states of the quantum field. 

There exist relevant antecedents to our analysis in the context of quantum field theory in a fixed curved spacetime, 
i.e., when the back reaction of the quantum fields on the spacetime geometry is not taken into account. The so- 
called Bunch-Davies vacuum [3l|, [s^l for fields in de Sitter is a state invariant under all the isometrics of de Sitter 
space, which is maximally symmetric. The renormalized expectation value of the stress tensor operator for that state 
is proportional (with a constant factor) to the metric and therefore its contribution to the semiclassical Einstein 
equation has the same form as a cosmological constant term. More importantly, it was shown in Ref. |35] that for 
fields with a wide range of mass and curvature-coupling parameters evolving in a given de Sitter spacetime, the 
expectation value of the stress tensor for any reasonable initial state tended at late times to the expectation value 
for the Bunch-Davies vacuum, where by reasonable states one means states with the same ultraviolet behavior as the 
Minkowski vacuum, i.e., with essentially no excitations at arbitrarily high frequencies (technically they arc known 
as Hadamard states [3l|, [S^l)- This result can be intuitively understood as follows: the exponential expansion will 
redshift any finite frequency excitations of the Bunch-Davies vacuum so that their contribution to the stress tensor 
will tend to zero at late times. 

The result described in the previous paragraph suggests that even when taking into account back-reaction effects, 
perturbations around de Sitter will be redshifted away and at late times the spacetime geometry will approach de Sitter 
space with an effective cosmological constant which includes the contribution from the expectation value of the stress 
tensor for the Bunch-Davies vacuum. However, in order to prove this expectation without any room for doubt, one 
needs to solve both the semiclassical Einstein equation and Klein-Gordon equation for the scalar field self-consistently. 
This is the main goal of this paper. We will consider a fairly general family of Gaussian initial states for the quantum 
field which are spatially homogenous and isotropic, and discuss under what conditions the trace of the stress-tensor 
expectation value exhibits unphysical divergences at the initial time after the standard renormalization procedure. 
We will also explain how to select appropriate states with regular initial behavior, and then solve the back-reaction 
equation explicitly. The standard renormalization procedure for the ultraviolet divergencies of the expectation value 
of the stress tensor requires the renormalization of the gravitational and cosmological constants as well as two new 
dimensionless parameters which are related to local geometric terms in the gravitational action which are quadratic 
in the curvature tensor. These new parameters should in principle be determined experimentally in order to eliminate 
the two-parameter ambiguity otherwise exhibited by the back-reaction equation. Nevertheless, for the particular case 
that we are considering (and at the order in the Planck length at which we are working) the results turn out to be 
independent of the particular value of these renormalization parameters. 

The plan of the paper is the following. In Sec.|lT]we introduce the particular model that we will be considering and 
describe the procedure that we will employ to generate the initial state of the quantum field. In Sec. IIIII we use the 
Closed Time Path (CTP) formalism, whose application to semiclassical gravity is briefly summarized in the Appendix, 
to derive the effective action which encodes the effect of the quantum fields on a spatially flat RW spacetime. We 
also derive the relation between the time-time component of the stress tensor and its trace, which will be later used 
to derive the back-reaction equation. In Sec. llVl the CTP effective action is used to obtain the quantum expectation 
value of the stress tensor operator, and the suitable initial conditions for the quantum state of the field and for the 
cosmological scale factor are discussed in detail. Finally, in Sec.|V]the semiclassical Friedmann equation that describes 
the back reaction of the quantum field on the scale factor, driven by a cosmological constant, is derived and solved 
perturbatively in powers of the Planck length over the Hubble radius. We find an explicit result valid for all times. 
It shows that there exists a self-consistent de Sitter solution with a slightly shifted cosmological constant due to the 
one-loop effects and that all the other solutions tend to this one at late times. 
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II. A QUANTUM FIELD IN A RW BACKGROUND AND ITS INITIAL STATE 



In this section we describe our model for the back reaction of quantum fields on a cosmological background. We 
assume a spatially homogeneous and isotropic cosmological model with flat spatial sections, described by the metric 

gt,v = a'^{v)Vt^-y, (1) 

where 77^^ is the n-dimensional Minkowski metric (we use arbitrary dimensions for the moment in order to perform 
dimensional regularization later on), which takes the form 77^^ = diag(— 1, 1, . . . , 1) when considering the usual inertial 
coordinates, and 0(7?) is the cosmological scale factor in terms of the conformal time r], which is related to the physical 
time thy adrj = dt. The classical action for a real massless scalar field ^{x) coupled to gravity is 

Smig^u, '^] = -\j d''^V~9 [g^'d^'^d.'^' + (Cc + y)R , (2) 

where the dimensionless parameters = {n~2)/[A{n~ 1)] (equal to 1/6 in four dimensions) and v give the coupling 
to the Ricci curvature scalar i?, given in this case by 

i? = 2(n-l)(-3+^-,). (3) 

where here and throughout the rest of the paper overdots denote derivatives with respect to the conformal time, i.e., 
'= d/drj. The minimal coupling case (no direct coupling to the curvature) corresponds to = — a massless scalar 
field with minimal coupling mimics the behavior of gravitons in the cosmological background, except for a factor of two 
corresponding to the graviton polarizations. When v = Q the classical action S-[a[g^v, is invariant under conformal 
transformations with g^y — + Vl^{x)g^i, and <I>(x) — > f2^^~"^/^(a;) $(a:); this is known as the conformal coupling case 
and it can be used to mimic the behavior of photons. 

Since the RW metric given by Eq. ^ is conformally flat, it is convenient to introduce the rescaled scalar field 
(j){x) = a^"~^^/^(?7) $(a;). The action in Eq. ([2]) then simplifies to 

5m[a>] =-\j d''x{v^''d^cj)d,<l> + iya^R^^) , (4) 

which is the action for a free scalar field (j){x) in Minkowski spacetime with a time-dependent quadratic coupling 
va^Rcj)^. Identifying the matter Lagrangian from S'm — J drjL^ and the momentum Tr(x) — 5L^/ 5(j){x) — 4>{x), 
one obtains the following Hamiltonian for the rescaled scalar field: 



TT^ + (vA^ + va^R(f^ 



H^^+Hm, (5) 



where in the second equality we have separated the Hamiltonian into that of a free massless field, i?in\ and an 
interaction Hamiltonian iJint, which is proportional to v. 

Let us now discuss the kind of initial quantum states of the field that we will be considering. We are interested in 
the evolution of the scale factor driven by a cosmological constant A plus the back-reaction effect due to the quantum 
scalar field, given some initial conditions for the scale factor and its derivative at some initial time r]i as well as the 
initial state of the quantum field at that time. On the other hand, we will use the evolution from 7/ — —00 to 7/ = iji 
as an auxiliary way to prepare the initial quantum state of the field. More precisely, the initial state of the field will 
be given by 

1*0 =?7(7,i,-(X3)|0,-(X3), (6) 

where |0, —00) is the usual Minkowski vacuum at 77 — > —00 and U is the time evolution operator associated with the 
Hamiltonian given by Eq. ([5]). U{r]i, —00) only depends on the scale factor before the initial time rji, which will be 
denoted from now on by a$ (77) and can be a fairly arbitrary regular function, subject only to condition 

lim a|(77)i?*(77) =0, (7) 

77 — > — QO 



where i?* (77) is the Ricci scalar corresponding to a>ii (77) , and to the requirement of a sufficiently smooth transition at 
77i to the scale factor at later times (the reason for this latter condition will be explained in detail in Sec. lIVp . Thus, 
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our initial state is a squeezed state-'^ that evolves (in the Schrodinger picture) from the Minkowski vacuum state. A 
particular case is the Bunch-Davies vacuum [s^ for zero mass and curvature coupling (^c + i/), which would follow 
from considering a scale factor a*(ry) that corresponds to a given de Sitter spacetime all the way from rj — —oo to 
ij = rji, [note that such a scale factor does satisfy condition ([7])]. Note that the state obtained from the construction 
described above and defined by Eq. ([6]) is the state of the rescaled field <j>{x). However, the state of the original field 
^(x) can be derived straightforwardly from it if one takes into account the simple relation between (l){x) and ^(x) 
involving the scale factor 0(77). 

Finally, the gravitational action for the scale factor can be written as 

Ss [«] = ^ drj a^\R- 2 A) + 5| [a] , (8) 

where the first term is just the Einstein- Hilbert term with k = SttG = Sir /nip (G is the gravitational coupling 
constant, nip is the Planck mass and we are using natural units with h = c = 1), V = J d"~^x is a spatial comoving 
volume factor that will drop in the final expressions, and S^[a\ accounts for the gravitational counterterms that will be 
specified later. At this point the parameters k and A should in principle be considered as bare parameters. However, 
for a massless field these parameters do not need to be renormalized when using dimensional regularization since the 
divergencies in that case only require counterterms which are quadratic in the curvature, as we will see below. 



III. THE EFFECTIVE ACTION FOR THE COSMOLOGICAL SCALE FACTOR 
A. The expectation value of the energy density and the trace of the stress tensor 



We can now follow the procedure outlined in the appendix to derive the semiclassical Einstein equation describing 
the back reaction of the scalar field on the spacetime geometry. The expectation value of the stress tensor for a given 
state of the field plays a key role in that equation. As long as one considers spatially homogenous and isotropic states 
of the quantum fields, it is consistent to assume that the metric in the semiclassical equation (jA.l[) takes the restricted 
form ([T]) throughout, so that there is only one dynamical variable 0(77) to be determined. 

Hence, we can concentrate on just one of the equations for the different components of Eq. (jA.ip . and in particular on 
the 00 component, which in this semiclassical cosmological context may be called the semiclassical Friedmann equation. 
The expectation value (Too)rcn will be taken in the state defined by Eq. ([6]). Both the classical stress tensor and 
its quantum expectation value can be obtained by functionally differentiating, respectively, the classical action Si-a 
and the influence action S\f with respect to the metric, according to Eqs. (|A.7[) and (jA.Sp . The influence action 
describes the effect of the quantum matter fields on the gravitational field and results from functionally integrating 
the quantum matter fields. Note, however, that since our metric has been assumed to have the from g^i, — a'^{ri)'r]^^, 
with the scale factor as the only independent kinematical degree of freedom, we can only functionally differentiate 
with respect to 0(77) and will just be able to obtain trace of the stress tensor, T^. This can be seen as follows. Since 
^gfj,L/ = 2arj^i,6a = 2a~^g^v5a, from Eq. (|A.7p we can write 



55„, = j d''x^6g^.= I d"xV^T^'''g^,^ = j dii j d 



(9) 



so that 5Sni/5a — Va"~^T^, where the spatial comoving volume V appears as a consequence of the spatial homogeneity 
of the stress tensor. This also means that by functional derivation with respect to the scale factor we obtain, from the 
renormalized influence action S*™ in Eq. (jA.Sp . the renormalized expectation value of the trace of the stress tensor 
operator: 



ss 



IF 



(10) 



where from now on we will drop the subscript "ren" in the expectation value to simplify the notation. 

In principle one could use either the 00 component of the semiclassical Einstein equation (the Friedmann equation) 
or the equation for the trace. There is, however, a subtle difference. For the sake of the argument let us consider the 



This means that we restrict our attention to a family of Gaussian pure states. 
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classical limit only; whereas the classical Friedmann equation is a first-order differential equation in time for the scale 
factor the equation for the trace is of second order in time. The solutions of the Friedmann equation automatically 
satisfy the trace equation, but the Friedmann equation also constrains the initial conditions for the scale factor and 
its time derivative. If one works only with the trace equation, this additional information is missed. (When quantum 
corrections are added, higher-order derivatives appear in both equations, but we will explain how to deal with them 
in Sec. M) 

In this paper we will work with the semiclassical Friedmann equation. Therefore, we need to calculate (Too). One 
possibility is to start with a metric of the form ds^ — —N'^{ri)dTf + a^{jj)5ijdx^dx^ with N{fi) and a{fj) independent,^ 
functionally differentiate with respect to N and a, and finally take N — a only after that. The functional derivative 
with respect to N gives (Too) and the Friedmann equation. Alternatively, one can make use of a useful relation 
between (Too) a-nd (T^) in a RW spacetime which is a consequence of the stress-tensor conservation law V^(T'''') — 0, 

and the fact that ^= d/drj is a conformal Killing field, i.e., 2V(a'fb) = ^g^v with A — 2d/a in our case [3ll,[3H|- These 
two equations lead to 

v^((fne.) = ^(r;'), (11) 

which can be integrated over the spacetime volume bounded by the space-like hypersurfaces corresponding to ry' = — oo 
and rj' = T] (for the construction introduced in Sec. |TT] to generate the initial quantum state of the field the stress 
tensor is also conserved from rj' — —oo till our initial time rj' = rji). Using Gauss's theorem we get a relationship 
between the integration of (Too) on the two hypersurfaces and the spacetime integral of (T^) which reads 

{fM)a--\l)^C-f dry'a"-i(r/)d(r;')(r; (,?')), (12) 

where C = (Tbo(— oo))a"~^(— oo) and we have divided by the spatial volume V, which appears due to the spatial 
homogeneity of the stress-tensor expectation value. Since the constant C is proportional to the expectation value of 
the energy density of a Minkowski vacuum, as follows from Eq. ([S]), it should vanish. Hence, from now on, we will 
take C = 0. 



B. The effective action 



In this subsection we compute the influence action needed to derive the expectation value of the trace of the stress 
tensor according to Eq. pil)) . Using the so-called closed time path (CTP) formalism, the influence action Syplg^.g^] 
for an arbitrary metric is defined in the appendix by Eqs. (IA.4[) or (jA.lOp for a general initial state of the field. 
Specializing Eq. (jA.lOp to the conformally flat metric g^i, = a?"q^i, with the definition of our initial state as given in 
Eq. ([6]), we get 

g»S:p[a+,a-] ^ / (0, -(X3 1 J7_ (-00, ??( ) | ^) (^| J7+ (r/f , -oo) |0, -(X)) , (13) 



where \ip) are the properly normalized field eigenstates, such that (p{x)\Lp{x')) = Lp{x)\ip{x')) , and the time evolution 
operator is 

lf±{r]f, -oo) =Texp^-i J dr]H^[a±,TT,(t)]^ . (14) 

Using the path integral representation for the time evolution operator in terms of the action ^ for the scalar field 
we have 

g.SiF[a+.a-] ^ f p<^+p^-g»(S„.[a+,0+]-S;;[a-,0-])^ (-^5^ 



^ Provided that N{'r}) is non-vanishing and differentiable, such a metric can always be rewritten as ds^ = a? {ri){—d'q'^ + Sijdx^ dx^ ) through 
a coordinate transformation involving a redefinition of the time i). 
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where (^^{rjf) = (f>~{rif). This expression corresponds to Eq. (|A.4[) for an initial Minkowski vacuum state. To enforce 
this state we must take the usual —ie prescription. Integrating by parts and taking into account that the action 
is quadratic in the field, we can write 



(det A) 



-1/2 



(16) 



where the matrix A is defined by = r/^'^9p9iy — i^(a+)^i?+ + ie, j4__ — — yrj^'^ di^id^, — v{a^y'R^ — ie), and 

A^ = A ^ = 0, and a Gaussian integration has been performed in the last equality. Introducing the inverse matrix 

G = A^^ we thus have 



5iF[a+,a ] - 



trlnG. 

2 



(17) 



The matrix G can be computed perturbatively. Following Refs. [13, [H,!!^ we define A = A'^ + V, where the matrix 

V includes the time-dependent interaction with V^^^^ — — i^(a+)^i?+, V = i'{a~)^R~ . Then up to second order in i^, 

G = G°{1 - VG° + VG°VG° + ...) where G° is the Minkowski 2x2 CTP propagator with G% = Ap, G° „ = -A^,, 

G*^ = — A+ and G'L^ = A^, and where Aj;- and A^i are, respectively, the Feynman and Dyson propagators and 

A^ are the Wightman functions: 



(27r)" p2 ^ie' 

Substituting into Eq. p7p we have (up to second order in v) 



(18) 



SiF[a+,a ] 



--trlnG" + -tr(F++A;^) - -tr(T/__Ac) 
-UTiV++AFV++AF) ~ -lT{y—ADV-^AD) 
+ ^tr(y++A+F__A-). 



(19) 



The first three terms do not contribute to the dynamical equations for a(r\), the first term is independent of a, and 
the second and third terms are tadpoles which are identically zero in dimensional regularization |;4ft] . so that there 
is no linear term in v in the effective action. The fourth and fifth terms involve the product of Fey nman and Dyson 
propagators and need regularization, whereas the last term is finite. Following closely Refs. [371 [sst we get, after 
expanding in powers of (n — 4), that the real part of S'lp in n dimensions is 



ReS'iF , a 




dT^dri'A ( - (?7) ) #(?7 - 77') S ( - (?7') 
a I \a 



+ 0{n-A), 



(20) 



where we have used the difference and semisum notations A(/) 
and the kernel H{rj — rj') is given by 



/+ - /- and S(/) = (/+ + f-)/2, respectively. 



2^ 



ZTT 1 

In \uj\ + ysign(-cj) - -(2 + ln47r - 7) 



(21) 



where 7 is the Euler-Mascheroni constant. Note that although the argument of the logarithm is not dimensionless, 
when combining the influence action with the counterterms in the bare gravitational action, the contribution involving 
the renormalization scale /i will finally render the argument of the logarithm dimensionless, as we will see below. There 
is also an imaginary part in Sip, but it does not contribute to the expectation value pH)) . and thus to the semiclassical 
equation for a, because it depends quadratically on the difference variable A(a/a). This means that when functionally 
deriving with respect to a"*" and then taking a+ = = a to get the expectation value, the imaginary contribution 
vanishes, as it should. The role of the imaginary part of S'lp is related to the so-called noise kernel, which accounts 
for the fluctuations of the stress tensor and allows to go beyond the semiclassical equations, which capture only the 
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averaged value of the stress tensor. The noise kernel plays a key role in stochastic gravity, see Refs. [4l|, |42, [43, |4j| 
for the general theory and Refs. [H, H^, IH, [13, El] for cosmological applications. 

As explained in the appendix, the dynamical equations for the gravitational field can be derived from the so-called 
CTP effective action, r[a+, a~], the two ingredients of which are the gravitational action [a], given by Eq. ([H]), and 
the influence action 5iF[a'^, a^]- Specializing Eq. (jA.5|) to the conformally flat metric in Eq. ([T]), the regularized CTP 
effective action becomes 



So; 



+ 5iF[a^, a ]. 



(22) 



The real part of the regularized influence action is given by Eq. (|20|) and diverges for n = 4. Thus, we need to add 
appropriate local covariant counterterms, which we denoted by 5*1 [a], to the bare gravitational action. For a massless 
field only counterterms quadratic in the curvature, as explicitly given by Eq. (|A.12[) . are needed. The integrand of the 
first term on the right-hand side of Eq. (|A.12[) . which is independent of i^^, is proportional to the square of the Weyl 
tensor in n = 4 dimensions. For a conformally flat metric like the metric ([1]) that we are considering here, the Weyl 
tensor vanishes and thus this term vanishes in four dimensions. However, it plays a crucial role in the trace anomaly 
[Slj . In fact, expanding in powers of (n — 4) we have 



a J \ a 



+ 0((n-4)2). 



(23) 



which is of order 0{n — 4) and therefore gives a finite contribution when multiplied by the divergent (n — 4)^^ factor. 
On the other hand, due to the 0{n — 4) dependence there will be no contribution proportional to the parameter a 
from Eq. (jA.12p . as expected since the tensor in Eq. (|A.1|) vanishes for a conformally flat metric. As for the 
second term on the right-hand side of Eq. (|A.12[) , we have 



,n—4 



+ 12{n~A) { 3 



ln(a /i) 



+ 0((n-4)2). (24) 



When multiplied by (n — 4)~^, the first term of this expansion cancels out the (n — 4)~^ divergence of Eq. (pOj) . 

Finally, we can add the regularized counterterms of Eq. (|A.12p with the particular values in Eqs. ((23l) and (jM)) to 
the Einstein-Hilbert action (including the cosmological constant) to obtain the total bare gravitational action S'g[a]. 
Together with the regularized influence action, whose real part is given by Eq. (|^D|) . it gives the regularized CTP 
effective action r[a"'", a~]. We can then take the limit n — > 4 to obtain the four-dimensional effective action in terms 
of the renormalized gravitational action and influence action: 



rcn + 



a ,a \. 

The result for the renormalized real part of the influence action in four dimensions is: 
ReS'[^"[a+,a-] I drj ' 



(25) 




A( (^] Ina -2 /d^'Ar(r,)]H(r,-,/;/2)Er(^') 



where we have incorporated the renormalization scale fj, — /iexp[(2 -f ln47r — 7)/2] in the new kernel 



—e ^^^^ " M In i,r^ + — sign(-tj, 

ZTT \ H Z 



(26) 



(27) 



From the action in Eq. (|26p one can obtain the renormalized expectation value of the trace of the stress tensor 
(remember that the imaginary part of the influence action plays no role in that). Similarly, by taking the functional 
derivative of the CTP effective action r[a"'",a~] with respect to and then equating = = a we obtain the 
trace of the semiclassical Einstein equation. Note that the CTP effective action, given by Eq. (|22|) or equivalently by 
Eq. ()25p . is renormalization-group invariant, i.e., it is independent of the renormalization group scale and so are the 
physical predictions that one can derive from it. The dependence on /x in ReS'[p"[a''", a~], which gives rise to a local 
term of the form (Qi^^/Stt^) In/Lt A((d/a)^), is exactly compensated by the the dependence on /i of the renormalized 
parameter (3 multiplying the B? term in the renormalized gravitational action. This can be traced back to the fact 
that the bare parameter /3b is independent of /x, as explained in the appendix. 
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We close this section by mentioning an alternative (but entirely equivalent) method of calculating the influence 
action provided in Ref. The approach, which is based on decomposing the field in spatial Fourier modes, 

computing the unitary evolution operator for each mode perturbatively in the interaction picture, and summing over 
all the modes at the end, can be useful when considering more general initial states at a finite initial time r]i which 
are not necessarily of the form given by Eq. ([6]). 



IV. THE EXPECTATION VALUE OF THE STRESS TENSOR 



A. The trace 



Functionally differentiating with respect to a+ the expression for the influence action given by Eq. (pS)) and using 
Eq. (fTO]) . we obtain the expectation value of the trace of the stress tensor: 
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(28) 



where t]) ~ drj' H {rj — ij' ; n)f{ri'). One can easily check that the usual result for the trace anomaly is obtained 
in the conformal limit h' ^ 0. Indeed, taking into account that in four dimensions 



dr/"^ 



and 

Eq. ((28l) can be rewritten as 

in 



dr]\a'' ~^ " 
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28807r' 



■UR- 
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28807r2 



{R^.p.R^'""' - Rf^.R^ + 0{y^), 



(29) 



(30) 



(31) 



which coincides with the trace anomaly [3l| for a massless conformal scalar field when v — 0. Note that the counterterm 
in Eq. ((23)) . when multiplied by the divergent factor (71 — 4)^^, plays a key role for this result. 

Eq. (|28p has a non-local term which includes the functional ^[a/a; 77) and its first and second derivatives. Let us 
examine this non-local part in some detail. The Fourier transform in Eq. P7)) can be computed [see Eq. (VII. 7. 18) in 
Ref. [H] to yield 



0{ri - rj 

v-v' 



(7 + ln^)(5(?/-?7'), 



(32) 



where Vf stands for Hadamard's finite part prescription, and 7 is the Euler-Mascheroni constant. This prescription 
means that 



k[/;77) 



lim 



drf 

00 v-v 



-/(77') + (lne-l-ln/H-7)/(77)| 



(33) 



Since the evolution from rj — — cx) to T/i can be regarded as an auxiliary way to generate the initial state of the quantum 
field, which is determined by the scale factor a(r]) for times 77 < rji and denoted earlier by aqi{ri), it is convenient to 
define ^(77) = {d^, /a^^,){r|). Therefore, for rj < rji we have K[d^ /axi,;rj) = K[v;ri), whereas for rj > rji the integral in 
Eq. (|33|) can be separated into two parts: 
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where 



— ; ?7 ) = — lim 
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K2 [v.Vl) , 



'—-[t]') + Ine + In^ + 7-77 

rj — rj' a a 



(34) 



(35) 
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which involves a time integration only after 7]i, and 

K2[v;r])^-[ , v{-q'), (36) 

which involves a time integration only before r]i. Following this notation one can also separate the trace in Eq. ([^5]) 
for rj > rji into two parts: (T^) = (r^)i + (T'^)2, where the first term involves Ki[d/a;ri) and the second one 
involves K2[v;ri). From Eqs. ([5]) and ([6]), and taking into account that in four dimensions a^R — 6a/a, we see that 
v{ri) completely determines the initial state of the quantum field. Thus, (2^)1 is state-independent, whereas (2^)2 
contains all the dependence on the initial state. In fact, the former will appear in the trace even if the initial state 
does not have the form given by Eq. 



B. The initial conditions 



We can now see that the state-independent part of the trace of the stress tensor, namely (T^Oii diverges in the 
limit rj — )■ i]^ , corresponding to the initial time. Let us consider the definition of Ki[d/a;r]) in Eq. (|35|) and Taylor 
expand d/a around rj' = r/ so that (a/a)(r/) = (a/a)(r/) + 0(77 — rj'); we then have 



Vi 



^ ^ d d d 

- iv') = — (^) In e + - iv) \n{^-v,)+0{v-m), (37) 
•q — T] a a a 



which implies that 



-■tA =--(?7)[ln(77-r;i)+ln/l + 7] + 0(7?-r;i), (38) 
a I a 



so that Ki\d/a]rj) is finite for all 77 > 771, but diverges when rj rf^ . Moreover, in addition to K,i[d / a; r]) , (T^')i also 
contains the first and second time derivatives of ki. Obviously upon derivation Eq. (|38p becomes more singular at 
77 — !■ r/j^, and in general the divergencies of the three singular terms will not cancel out. Therefore, — > 00 in the 
limit 77 — !■ 77j^. 

Thus, (T^) will diverge at the initial time unless we choose an initial state that cancels the previous divergencies. 
This partly motivates the class of initial states we consider in this paper. Taylor expanding v around rj' = rji, and 
following a procedure analogous to the one that led to Eq. ([HT]) [notice that condition ([7]) ensures that K2[v;ii) is well 
behaved in its lower limit], one can check that 

K2[v; rj) — v{r]i) hi{r] — rji) + finite terms, (39) 

whereby "finite terms" we mean terms which are finite in the limit 7; 771. Therefore, k [a/a; 77) — Ki[d/a;ri) + K2[v;ri) 
is finite at rji provided that f(77i) — lim^^j^+ (a/a). If this condition is satisfied, it is easy to see from Eqs. ([55]) and 

()36|) that ^^[a/a; 77) = K[^(a/a); 77). Thus, by iterating the same argument we find that the conditions that an 

initial state of the class ^ must satisfy in order to avoid initial time divergencies of (T^) are 

-y(77i) = lim ( - ) , 77(771) lim ( - ) , i;(77i) = lim ( - ) . (40) 



ya/ drj \aj di]'^ \a 



In other words, the cosmological scale factors before and after the initial time must be matched with continuity up to 
the fourth derivative. Eqs. (|40p have been presented as conditions on the preparation of the initial state given a scale 
factor a{ri) for rj > rji. Alternatively, one can regard these equations as initial conditions on the scale factor given 
some initial state defined by Eq. Henceforth we will assume that these conditions are satisfied. 



C. The energy density 



Once the expectation value of the trace of the stress tensor is known, the 00 component of the stress tensor may 
be obtained from Eq. The local part of aa'^(T^) is a derivative, so it gives rise to local terms in a^(roo). The 
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result, for ry > 771, is 



Hfoo) = 
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28807r2 
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(41) 



where T, which includes the non-local part of this expectation value and comes from integrating the last two terms 
in Eq. is defined by 



T[a, v;ri) — / drj' a 



drj' 
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(42) 



Making use of the "smoothness" of the scale factor at rji assumed above, and following the previous notation for the 
separation of the non-local terms, we can again write 



T[a,v] = Ti[a] +T2[a,v], 



(43) 



where we split r[a,f] into two parts: one independent of v and a second one which contains the entire dependence 
on V. This is achieved by splitting the integral in Eq. (|42l) as well as the functional k. The first part is given by 



Ti[a;ry) 



Kl 



a , 
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+ 2 - (r/) 
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a ' 



Kl 



drj'' 



- ;v 

a ' 



(44) 



This is obtained from Eq. (|42)) by explicitly applying the second-order derivative d^/drj'^ and taking into account that 
when the conditions in Eqs. (|40p hold, the derivatives acting on the functional k can be taken inside and applied to 
the argument. All this is done before splitting the integral and the functional n. Similarly, the second part in Eq. 
can be written as 



T2[a,v;ri) 



dv' -]{rf)\2 



'dri' 



K2 [v;ri 



+ - ('7i) i^[v;m) 



a ' 



(V) K2[v;'n') - K2[v;v' 
dri'v{ri')K[v; 77'). 



(45) 



The different form of the last three terms, which correspond to the integral from —00 to ?7i, is because we made use 
of the following equivalent expression for the last two terms in Eq. (pS)) times {—ho?): 



drj' 



1 



drj' \ drj' 
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a d 
a drj' 



(46) 



In addition to being convenient later on, we did this so that it became manifest that T2[a,v] depends on only 
through V. From Eqs. (P^ - P5|) it is clear that all the dependence of the energy density on the initial state of the 
quantum field is included in T2[a,v]. We should, however, remember that conditions have been used so that 
the divergencies from the evolution after rji and those from the initial state would cancel out. Consequently, the 
integrand on the right-hand side of Eq. and the integrand of the first integral in Eq. (|l5)) only exhibit logarithmic 
divergences when 77' — > rji which would cancel out when adding Ti and T2. In fact, they give a finite contribution even 
separately since, being logarithmic, they are finite upon integration. 

The energy density that we have just computed is valid for any scale factor at rj > rji and an initial state corre- 
sponding to any regular function v{rj) which satisfies conditions ([7]) and ([3D]). It can be included in the semiclassical 
Friedmann equation in the presence of any other classical source. In this paper we concentrate on the case where the 
classical source is a cosmological constant. 



V. THE SEMICLASSICAL FRIEDMANN EQUATION 
A. The equation 



Next, we proceed to analyze the evolution of the scale factor after the initial time rji, when it is driven by a 
cosmological constant and a free massless scalar field non-conformally coupled to the curvature. We will use the 00 
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component of the semiclassical Einstein equation (|A.1|) . namely, Gqq = —Agoo — (3Boo + K(T'oo)rcn- Note that we have 
already taken into account that A^,y = in our conformally fiat spacetime. On tlie otlicr hand, evaluating B^^, as 
given by Eq. (|A.3p . for the RW metric ((D), we get 



72 
7^ 



(47) 



which can also be obtained by using Eq. (j24p for n ~ -i and noticing that the term proportional to A (a/a) gives rise 



to the first term on the right-hand side of Eq. (|¥T|) . Taking into account that Goo = 3(d/a) 
using Eq. (|4ip we obtain the following expression for the semiclassical Friedmann equation: 



for the metric ([T]) and 



= H 



2„4 




drj 



3607r 



18i/^ 



In a 



dr/ 



- T[a,v] 



(48) 



where /? — 5767r/3 — l/(607r), and we have introduced the Hubble constant H = (A/3)^/^ and the Planck length 
Ip = VG — (K/87r)^/^. These two constants introduce two different time-scales into the problem, the Hubble time 
and the Planck time tp = ^p, which we will assume to be well separated, namely, ^ Ip. It should also 

be emphasized that the semiclassical Friedmann equation (|48l) is invariant under conformal-time translations in the 
following sense. We already know that v{ri), with domain (—oo^rji], characterizes the initial state of the quantum 
field, |5'(?7i)), when it is of the form given by Eq. The time translation of f (77), given by v^{ri) = v{rj + A), is 
defined in the domain (—00, 771 — A] and corresponds to an initial state |5'A('7i — A)) = |4'(?7i)). It is then easy to see 
from Eq. (|48p that if a(ry) is a solution for some initial state \^{ri\)) characterized by w(r/), then a{ri + A) is a solution 
for the initial state |^A('7i — A)), characterized by va{v)- 

Our semiclassical Friedmann equation (|48p . which is a non- linear third-order integro-differential equation, looks 
like a typical back-reaction equation. Because of the higher-order time derivatives, those equations exhibit extra 
degrees of freedom which usually translate into unphysical runaway-type solutions [50|. [Indeed, if we had proceeded 
analogously to what was done in Ref. [s^ for QED, the Friedmann equation and the equation for the trace would fix 
the third and fourth derivatives of the scale factor at the initial time, given some freely specified initial values for a, d 
and a. This is in contrast to the classical case, where only a (or alternatively d) can be specified independently at the 
initial time.] To get rid of the unphysical solutions several methods have been proposed. In some methods one looks 
only for analytic solutions, in a suitable perturbative parameter, in order to select physical solutions only [s^. Issl. [56| . 
In some other methods the equation itself is changed in order to get rid of the higher-order derivatives; this is the 
case of the so-called order reduction method [13] ■ See Ref. [s^l for a review of the advantages and the shortcomings 
of the different methods in the context of semiclassical gravity. 

Given how we have proceeded in this paper, it is more natural in our case to interpret the point discussed in the 
previous paragraph as follows. The solutions of our semiclassical Friedmann equation (|48p must satisfy the the three 
initial conditions in Eqs. (j40p . This means that given some initial state, completely characterized by the function v{ri), 
the semiclassical Friedmann equation together with the equation for the trace fix 0(771) and 0(771) since the second, 
third and fourth derivatives at the initial time are fixed by Eqs. (j40p . Thus, there is in general a unique solution 
compatible with a given initial state. Certain initial states, however, give rise to solutions with characteristic time- 
scales of the order of the Planck time (corresponding to exponential growth or oscillatory behavior) , which lie beyond 
the regime of validity of the low-energy effective field theory approach that we have implicitly been using. In those 
circumstances the higher-order corrections involving terms with positive powers of the curvature and suppressed by 
the corresponding powers of the Planck mass can no longer be neglected and the low-energy expansion breaks down 
(since there are in principle an infinite number of such terms). We will consider only situations with no Planckian 
features where higher-order corrections are negligible. Hence, in the spirit of the effective field theory approach (valid 



for IpH <C 1), we will look for perturbative solutions'' in powers of (IpH) as 

a(77) = ao(?7) + {lpHfai{r^) + 0{{lpHf 



(49) 



A perturbative expansion may sometimes miss the right long-time behavior of the semiclassical solution. This can happen when the 
effect of the quantum corrections, although locally small, builds up over long times giving rise to substantial deviations from the 
classical solution. An example of such a situation is the evolution of a black hole spacetime when the back reaction of the emitted 
Hawking radiation is taken into account. One possibility in those cases is to modify the back-reaction equation using an order-reduction 
procedure and then solve the resulting equation non-perturbatively [58| . In the cosmological case considered here one can argue that such 
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where gq (77) satisfies the classical Friedmann equation. From now on we will use the subscript to indicate the classical 
unperturbed values. Without loss of generality we can focus on the solutions airj) with 09(77) — —l/Hrj defined for 
ryi < ?7 < since any other well-behaved solution not involving Planckian scales is connected to one of these by a time 
translation. The function characterizing the initial state is given by v{ri) = vo{ri) + 0((^pi7)^), and the conditions in 
Eqs. (go]) imply 



2 4 

-2' ^o(?7i) = --3 



vt v\ 

Substituting the expression (j49p into the semiclassical Friedmann equation 
order perturbation 01(77): 

TT 



12 



(50) 

we obtain an equation for the first 



1 1 



2 1 

if 1720^772 



+r] T[ao,Vo] 



(51) 



which is a first-order differential equation. Notice that there is no dependence on the free parameter f3. This is not 
surprising if one takes into account that /3 is the coefficient of a term proportional to OR in the trace [see Eqs. ((28)) 
and (123)] and that for the classical background Rq = 4A is a constant (where Rq is the Ricci scalar for the classical 
background) . 

Let us now compute the non-local term r[ao,wo] in the previous equation. After substitution of oq into Eq. (I44p . 
the state-independent contribution of this term reads 

Ti[ao;77) = --^ + ^!-^ + ^ - (52) 
277* 77^77^ 77^*77 2rif 

which is independent of the arbitrary renormalization scale /i. The latter is again because when substituting the local 
part of Ki[/;77') into Eq. (|44p . one obtains an integrand proportional to DRq = 0; this point could also have been 
anticipated from the fact that the term proportional to In /i should have the same form as the term proportional to 
/3, which has been found to vanish above."* 

Substituting ao and 7;o into Eq. (I45p . we get the following result for the state-dependent part of T[ao,vo]r]): 



rj — T] 



rj' 77' rj' 



1 r N 1 r- ^ 

H — ^K[vo;'ni) H K[vo;r]i) + 



dj/vQ{f]')K[vQ;'q') 



(53) 



The conditions ([50]) guarantee that the dependence on fi of the last three terms cancels out. The integral over rj' in 
the first term can be easily computed. Finally, adding up the two contributions Ti and T2 we get 



T[ao,vo;r]) 
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„2 



where 



and 



drj' In 
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'7i 



dr]'vo{'i]')K[vo;'i]') 



(54) 



(55) 



(56) 



We have already pointed out that the first order equation (j5ip includes neither the free parameter P nor the 
arbitrary renormalization mass scale /z. Moreover, it does not depend on a either since = for a conformally 
flat spacetime. This means that the semiclassical Friedmann equation (j48p is fully predictive to this first perturbative 
order for the particular case of a de Sitter background that we are considering. 



an accumulation effect will not be present. That is because the classical cosmological constant implies a monotonous growing behavior 
for the unperturbed solution and the locally small effect due to the vacuum polarization of the quantum fields (much smaller than the 
classical cosmological constant) generates a perturbation which is always small compared to the unperturbed solution. Moreover, one 
can explicitly check a posteriori that the deviation from the classical solution does not become significant at late times. 
* Note that the results obtained in Refs. [26l . |27| in a similar context exhibit a non-vanishing dependence on /i. This should not be the 
case for the reasons given here. 
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B. The solution 



The equation (|51[) for the function 01(77) is a first order, linear differential equation, which can be easily solved 
in an explicit form. The general solution of a differential equation of the type y'{x) = —{2/x)y + f{x) is y{x) — 
k/x^ + (l/x^) / dxx'^f{x), where k is an arbitrary integration constant. Therefore, the solution of Eq. ([5T|l is: 



( 1 



27r 
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TTiJ 
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r*(77) 



where 



with 



2 2 1 

-;3^o(??') - —2^o{v') - —Mv') 
rj' rj' rj' 



- 77') 
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2 ^7 ^' 



1-^ 
rj' 



(57) 



(58) 



(59) 



The condition ([7]), which wo(?7) must also satisfy, ensures that the integral in Eq. ([55)1 is finite in its lower limit. As 
for the upper limit, all the possible divergencies in the integrand cancel out (moreover, even if they did not cancel 
out, they would give a finite contribution when integrated since they are logarithmic). Thus, rii,(77) is finite and, 
furthermore, it vanishes at future infinity, i.e., when 77 — > 0^. 

The solutions associated with the different possible values of the arbitrary integration constant k are related to 
each other by time translations: 0(77 + (/piJ)^A) — —1/Hr]+ (/piJ)-^ [01(77) + /S./Hrf] + 0((/piJ)**) . For the particular 
case fc = 0, and provided that H^^ 3> ?p, the first-order term in the expansion (j49p is much smaller than the 
zeroth-order term for all times 771 < 77 < 0, and one expects that the higher-order corrections can be neglected [since 
they are suppressed by an additional power of [IpHY] and the solution 0(77) can be approximated by the first-order 
perturbative expansion. This approximation is no longer suitable for solutions (initial states) associated with other 
values of fc, but we already know that they are just time translations of the solutions with k = Q. Therefore, we can 
conclude that for reasonable states, namely, those which do not involve Planckian scales (and up to time translations) 
the scale factor after the initial time rji is given by 



1 



,37^ 



B 



* 3 



where the modified Hubble constant is 



H = H 



1 



( 1 



V7207r 



r*(?7) 



27r 



0((/pi/)4) 



(60) 



(61) 



The result in Eq. ([50)1 is valid for all times (after the initial time). It is the sum of a corrected de Sitter solution 
plus terms that vanish at future infinity 77 — )■ 0^. The corrected de Sitter term dominates at late times, regardless of 
the initial state. Therefore, we conclude that de Sitter spacetime is stable under spatially-isotropic perturbations in 
semiclassical gravity. 

The semiclassical shift of the Hubble constant does not have a definite sign, it depends on the coupling between 
the quantum field and the curvature. For instance, while the shift is positive for the conformal coupling case (1/ = 0), 
when V = —1/6, which is expected to mimic the effect of gravitons, this shift is negative, implying a small time- 
independent screening of the cosmological constant. The deviation from the conformal coupling, i>, is what makes 
the Hamiltonian ([5]) time dependent through the term involving the scale factor. This is what allows us to prepare 
a wide range of initial states by evolution of the m- vacuum with different scale factors 0^(77). In fact, the effect of 
this time dependence of the Hamiltonian can be interpreted as particle creation. Since the term proportional to v'^ 
in Eq. (pT|) is negative, one can say that at late times the effect of the created particles is to slow down the de Sitter 
expansion by a small amount, whereas the other vacuum polarization term, already present in the conformal case, 
has the opposite effect. 

A particularly interesting initial state is the Bunch-Davies vacuum [^'60(771)), which is obtained by evolution of the 
m-vacuum from —oo to r]i according to the scale factor of de Sitter spacetime or, in other words, with ^0(77) — 2/77^. 
In this case, one can check that A^ = B$ = r^(7;) = 0, so the solution ([60]) is just the de Sitter solution with a 
semiclassically modified Hubble parameter. Therefore, the Bunch-Davies vacuum together with de Sitter spacetime 
with the semiclassically modified Hubble constant constitute a self-consistent solution of the semiclassical Friedmann 
equation, to which other perturbed solutions tend at late times. 
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VI. DISCUSSION 

In this paper we have computed the one-loop vacuum polarization for massless non-conformal fields in a general 
spatially-flat RW background. This has then been applied to studying the evolution of spatially-isotropic perturbations 
around de Sitter spacetime when the back reaction due to a massless non-conformal field is self-consistently included, 
which corresponds to solving the equations of semiclassical gravity for this case. There is a self-consistent solution, 
associated to the Bunch-Davies vacuum for the quantum fields, with an effective cosmological constant slightly shifted 
from its classical value due to the vacuum polarization effects. Furthermore, we have found that this solution is stable 
under spatially-isotropic perturbations since the perturbed solutions tend to it at late times. It should be stressed 
that our results are independent of the particular value of the renormalization parameters a and /3, and therefore fully 
predictive (at the first order in l"^ at which we are working) . It should also be pointed out that our results seem to be 
at variance with those of Refs. (26l. [27l|. In this respect, it is important to keep in mind that one should consider the 
stability of the self-consistent solution with the shifted effective cosmological constant rather than that of the classical 
solution obtained in the absence of vacuum polarization effects. Furthermore, when comparing the perturbed solution 
to the self-consistent de Sitter "attractor" it is important to take properly into account any relative conformal time 
translation since that can give rise to a spurious growth in time of their ratio. This can be clearly illustrated by 
comparing two copies of the same self-consistent de Sitter solution with a relative conformal time translation. 

We were able to obtain explicit analytic results by using two approximations. First, we considered a perturba- 
tive expansion in the parameter i/, which characterizes the deviation of the curvature coupling parameter from the 
conformal case, and truncated the expansion to quadratic order. Hence, our effective action is exact through order 

. Second, we introduced a perturbative expansion in powers of {IpH)'^. Its purpose was to obtain a fairly accu- 
rate description for phenomena involving length-scales much larger than the Planck length while discarding spurious 
solutions involving Planckian scales, where the effective field theory approach that we have been using breaks down. 
Truncating such a perturbative expansion for the solutions (rather than doing so at the level of the equation of motion 
and then solving it exactly) can sometimes miss the right long-time behavior. However, as we discussed in Sec. [Vj 
this should not be the case for the situation considered here. 

There are a number of natural extensions or generalizations of our work in this paper. First, one could consider the 
other possibility for having weakly non-conformal fields, namely, the case of fields with conformal coupling but with 
a small non- vanishing mass such that ^ H^. In that case one should be able to proceed analogously to what we 
did here by computing the effective action perturbatively in rri^ through order m'* and then solving the equation of 
motion through order {IpH)"^. Second, the case of strongly non-conformal fields with a large mass ^ could be 
explicitly calculated using an adiabatic (or WKB) expansion [s^. The effective action in that case could be written 
as a local expansion of positive powers of curvature invariants suppressed by the corresponding power of A/^, a form 
which can be anticipated from local effective field theory arguments based on power counting and taking into account 
the relevant symmetries (diffeomorphism invariance in this case). On the other hand, the case of strong non-conformal 
coupling to the curvature could be treated by introducing a field-dependent conformal transformation relating the 
original Jordan frame to the Einstein frame, ^ where the curvature scalar does not couple to the scalar field {61l]. 

Third, even though the construction based on Eq. ^ can generate a fairly wide family of squeezed Gaussian states 
by considering a sufficiently general form of the auxiliary scale factor 0^(77), other approaches are needed in order 
to deal with all possible Gaussian states or even non-Gaussian ones. One possibility is to make use of the method 
developed in Ref. ^62] for a fixed background spacetime, which is based on the construction of fourth-order adiabatic 
vacuum states (one can show that the class of states that we have considered here are compatible with their approach 
and encompassed by the general class of states it can deal with). However, even if we restrict ourselves to the subclass 
of initial states generated by Eq. ^ , the procedure can be straightforwardly generalized to states with non- vanishing 
expectation values of the field or its canonically conjugate momentum. Since we are considering Gaussian states, this 
can be done by decomposing the field as a sum of a classical part which characterizes the evolution of the expectation 
value and satisfies the classical equation of motion plus a field with vanishing expectation value whose initial state 
can be generated by Eq. ([6]). This implies that when solving the back-reaction equation in powers of (IpH)^ one first 
needs to find the self-consistent solution for the background scale factor and the classical configuration of the field 
(corresponding to the evolution of its expectation value), and then solve for the perturbation of the scale factor due 
to the vacuum polarization effect of the quantum fluctuations evolving on that background geometry. 



^ Any physical predictions derived in the Jordan and Einstein frames should be equivalent, at least at the classical level. In the quantum 
mechanical case one still expects such an equivalence for small perturbations of the metric and the scalar field around their mean values 
[60l . In that case, however, it may be necessary to treat the perturbations of the metric and the scalar field on an equal footing, and 
the metric perturbations may need to be quantized. 
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Throughout the paper we have considered spatially- isotropic perturbations, i.e., RW geometries and quantum states 
compatible with their symmetries (spatial homogeneity and isotropy). It would be interesting to study the stability 
of de Sitter spacetime with respect to general inhomogeneous and anisotropic perturbations. In that case one has 
non-trivial results even for a massless and conformally-coupled field. The underlying reason is that when considering 
inhomogeneous and anisotropic metric perturbations, the perturbed geometry is no longer conformally flat. Whereas 
the stability of Minkowski spacetime with respect to general linear perturbations has been studied for arbitrary masses 
and curvature couplings (see Refs. [H, and references therein), such results do not exist for RW backgrounds. In 
order to analyze the dynamics of inhomogeneous and anisotropic perturbations around a spatially-flat RW background 
one can make use of the effective action and the semiclassical Einstein equation for general linear perturbations around 
a spatially-flat RW background obtained in Ref. [13] for a massless and conformally coupled field (see also Ref. 
for a more compact form). Unfortunately this linearized semiclassical equation is a cornplicated integro-differential 
equation and its solutions have not been studied in detail so far®; see, however, Ref. [64| for recent work in this 
direction. 

We have studied the back reaction of the quantum fields on the dynamics of the spacetime geometry within 
the framework of semiclassical gravity, which can be understood as a mean field approximation where the mean 
gravitational field is described by a classical metric whereas its quantum fluctuations are not considered. In order 
to study the quantum fluctuations of the gravitational field one can consider the metric perturbations around a 
background geometry corresponding to the semiclassical gravity solution and quantize them within a low-energy 
effective field theory approach to quantum gravity ,28, 29, 65., .66J. So far this approach has been mostly applied 
to weak-field problems [30|, but it seems particularly interesting to extend its application to str ong -field situations 
involving black holes and cosmological spacetimes (63, ■ The stochastic gravity formalism [i^ |44| | can be a useful 
tool in this respect since one can prove its equivalence to a quantum treatment of the metric perturbations if graviton 
loops are neglected, which can be formally justified in a large N expansion for a large number of matter fields [69j . 
A central object in this formalism is the symmetrized connected two-point function of the stress tensor operator 
for the quantum matter fields, which determines the metric fluctuations induced by the quantum fluctuations of 
the matter fields. Such an object has been computed for a massless minimally coupled field evolving in a dc Sitter 
background spacetime and the fluctuations of the stress tensor were found to be comparable to its expectation value 
[70| . Therefore, studying in detail the quantum fluctuations of the metric in this context constitutes a natural 
extension of our work worth pursuing. The results obtained here would still be relevant in that case because they 
provide the right background around which the metric should be perturbed and quantized. 

We close this section with a brief discussion of the relationship between our results and the linearization instability 
for metric perturbations around de Sitter spacetime coupled to a scalar field found in Ref. [7l| , where it was concluded 
that it is only consistent to consider de Sitter invariant states for the quantum field. This conclusion does not 
directly affect our analysis because we did not consider fluctuations of the metric and studied only the dynamics 
of the mean geometry, which couples to the expectation value of the stress tensor operator of the matter field. 
The expectation value of the stress tensor for the class of states that we have considered in this paper, which are 
spatially homogenous and isotropic, automatically satisfies the linearization stability constraint given by Eq. (44) in 
Ref. [71i |. It is when considering the quantum fluctuations of the metric that the linearization stability condition 
imposes additional restrictions on the state of the matter field because in that case the condition must be imposed on 
the n-point correlation functions of the stress tensor as well, and this implies that the state of the field must be de 
Sitter invariant. 
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For perturbations around Minkowski it is relatively easy to solve the integro-differential equation because it can be transformed into 
a purely algebraic equation by Fourier transforming with respect to not only the spatial coordinates but also the time coordinate. By 
contrast, that is not possible in the RW case due to the time dependence of the scale factor, and the non-locality in time cannot be 
eliminated. 



16 



APPENDIX: SEMICLASSICAL EINSTEIN EQUATION 



A 



ah 



In this appendix we briefly review the scmiclassical Einstein equation and its derivation by functional techniques. 
When neglecting the effects of graviton loops, the back reaction of quantum matter fields on the mean gravitational 
field is described by the scmiclassical Einstein equation, which can be written as 

Gab [g] + Apafc + aAab [g] + (3Bab [.?] = K{fab [.g])ron, (A.l) 

where Gab is the Einstein tensor associated with some globally hyperbolic spacetime with metric gab, {Tab [5]) rcn is 
the suitably renormalized expectation value of the stress-tensor operator corresponding to the scalar field operator 
(f> [g], and a, f3, A and k are renormalized parameters evaluated at the same renormalization scale as {Tab[g])icn- The 
tensors Aab and Bab are obtained by functionally differentiating with respect to the metric terms corresponding to 
the Lagrangian densities G'^'^^^Gabcd and in the gravitational action, 

^/ I d^x^gGcdefG^"^^, (A.2) 
V -5 ogab J 

B"'' = I d^x^R\ (A.3) 

V-5 ogab J 

where Gabcd is the Weyl tensor and R the Ricci curvature scalar. These Lagrangian densities CabcdG"^"'^ and R^ 
are related to the counterterms introduced in the bare gravitational action needed to renormalize the ultraviolet 
divergencies arising in the expectation value of the stress tensor. Note that from their definitions the tensors (|A.2|) - 
(IA.3|) are divergenceless: VAat = = VBab- 

Let us see how the scmiclassical Einstein equation can be derived by functional methods. To compute the expectation 
value on th e rig ht-hand side of Eq. (|A.ip we can use the closed time path (CTP) or "in-in" functional formalism 
[37I [67l [73, FtsTFtI. FtsI . F76| . Let us foliate the assumed globally hyperbolic spacetime with t = const, spacelike 
hypersurfaces St, and denote the initial and final times by ti and it, respectively. In the CTP formalism we introduce 
two copies of the metric and the field, (5^5, 5"^) and which will coincide at the final time: ff^;,(if) = gab^U) 

and 0+(if) = (j)~{tf). Let pi[(j)'^ (ti) , (ti)] be the matrix element of the density operator describing the initial state 
of the scalar field. The Feynman- Vernon influence action [73,[zH|, Swlg^, g~], which describes the effect of the matter 
field on the gravitational fleld, is deflned as the following path integral over two copies of the scalar fleld: 



where S',n[(7, (/i] is the action for the scalar field in the spacetime described by the metric gab- Neglecting graviton 
loops, the CTP effective action for the gravitational field is then 

T[g+,g-] - S,[g+] - S,[g-] + SMg+ , gl, (A.5) 

where Sg[g^] is the bare gravitational action. r[(7"'", g~] is the effective action for the mean gravitational field coupled 
to the quantum scalar fleld. SiF[g'^ , g~] has ultraviolet divergencies which can be renormalized by using a suitable 
regularization procedure and by adding the aforementioned counterterms to the bare gravitational action Sg[g^]. 
More specifically, one starts with a regularized gravitational action in Eq. (jA.SP which includes the bare parameters 
Kb, ^b,C(Bt Pb] at the end of the calculation one takes the regularization parameter to its physical value and the 
divergencies are absorbed into the bare parameters which acquire their dressed physical values. The renormalized 
effective action can then be written as 

T[g+, g-] - ^-"[5+] - Sl"^[g-] + ^[^[.9+, 5^], (A.6) 

where the superscript means that these terms have been already renormalized and are finite. 
Since the classical stress tensor of the matter field is defined as 

T-'' = ^^, (A.7) 
V-g Sgab 

one can see from the definition of the influence action in Eq. (|A.4[) that the expectation value of the stress tensor in 
the given quantum state of the fleld is given by 



59ab 



(A., 
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where the renormaUzed value of the influence action has been used. Finally, the semiclassical Einstein equations (|A.1[) 
can be derived by functional derivation of the effective gravitational action, r[g~^,g~], with respect to and then 
taking = g-f^ = gab'- 



Jab 



0. (A.9) 



Notice that doubling the number of fields, the "plus" field, which evolves forward in time, and "minus" field, which 
evolve backwards in time, is what allows us to obtain an expectation value from the above functional derivative in 
the CTP formalism, rather than a transition element as in the ordinary "in-out" effective action method. This can 
be clearly seen from the following alternative representation of the influence action: 



e 



iSiF[g+,s ] — ^j. 



Ug^{tt,U)piUl_{tt,U) ^('^i\Ug-{tutf)Ug+{tl,ti)\^i), (A.IO) 



where the last equality holds for a pure initial state of the field pi — and Ug± is the unitary time evolution 

operator for a field (j) propagating in a spacetime with metric g^^: 



tt 

± 



f/g±(tf,ti) =Texp -^ / dtH,^[g^,(t>,iT]] , (A.U) 



where T denotes time ordering and H^[g, (j3,TT] is the Hamiltonian of the scalar field obtained by considering the 
spacetime foliation {St} and introducing the corresponding 3 + 1 decomposition. 

In this paper we use dimensional regularization, this means that the regularization parameter is n — 4, where n is 
the number of spacetime dimensions. In that case the square of the Weyl tensor in the Lagrangian density of Eq. (jA.2[) 
must be substituted by | {Rabcd.R^'^ '^'^ ~ RabR'^^) , where Rabcd and Rab are the Riemann and Ricci tensors. Such a 
Lagrangian density reduces to the square of the Weyl tensor in n = 4 when the Gauss-Bonnet theorem is taken into 
account. 

We also consider a massless scalar field throughout. Using dimensional regularization the counterterms for a 
massless scalar field can be read from the divergent parts in the effective action computed using the DeWitt-Schwinger 
expansion of the Fcynman propagator in the coincidence limit; see Eq. (6.44) in Ref. [3l|. It should be remarked 
that the DeWitt-Schwinger expansion is defined for massive fields and it is ill defined in the massless limit. However, 
it can still be used to extract the divergent terms. There are three divergent terms in this expansion, the first and 
second terms have coefficients related to a constant and to the Ricci scalar, respectively. These two terms vanish 
in the massless limit. The third term has a coefficient which is quadratic in the curvature, and in this case it is 
convenient to keep the mass as an infrared regulator which can be removed at the end of the calculation by redefining 
the renormalization scale. This term leads to the following counterterm for the gravitational action: 

S|[g;^^] = {^(^') + 2880^(1 ~ 4) ) j d-x^g{RabcdR-'^' - RabR''') 



where a(/i) and are the renormalized dimensionless parameters which appear in Eq. (|A.1[) . and /i is a parameter 
with dimensions of mass which ensures that the action has the correct dimensions even when n 7^ 4 and plays the 
role of the renormalization scale in dimensional regularization. The bare parameters as and /3b, which correspond to 
the whole factor multiplying the first and second integrals respectively, should be independent of the renormalization 
scale PL. Therefore, under a change of the renormalization scale p p' [and neglecting terms with positive powers of 
(n — 4), which vanish when n — > 4] the renormalized parameters change as follows: a{p') — a{p) — (28807r^)~^ \n{p' / p) 
and I3{p') = Pip) - (j/V327r2) \n{p'/p). 

Note that as seen from Eq. the R^ counterterm in the action can be easily inferred from the (n— 4)~^ divergence 
in Eq. (|20[) . but it is not so obvious how to come up with the other counterterm (with the Riemann square) because 
no other divergence is present. That is because the first term in Eq. (|19|) . which would be the entire contribution 
for a conformal field, is finite in the particular case of a conformally flat metric. However, a slight departure from 
conformal flatness {e.g. by considering inhomogeneous perturbations around the flat space geometry) renders the first 
term in Eq. (fT9| divergent and such a divergence is accounted for by the above general DeWitt-Schwinger expansion 
[7^[80|. Hence, even though the first term in Eq. would still be independent of a when regulated with dimensional 
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regularization, the counterterms that cancel the divergencies in the limit n —^ 4 would bring the a dependence into 
the problem. 



[1] A. H. Guth, Phys. Rev. D 23, 347 (1981). 
[2] A. D. Linde, Phys. Lett. B 108, 389 (1982). 

[3] A. Albrecht and P. J. Steinhardt, Phys. Rev. Lett. 48, 1220 (1982). 
[4] A. D. Linde, Phys. Lett. B 129, 177 (1983). 

[5] A. D. Linde, Particle physics and inflationary cosmology (Harwood Academic, Amsterdam, 1990). 

[6] G. F. Smoot et al., Astrophys. J. 396, LI (1992). 

[7] C. L. Bennett et al., Astrophys. J. Suppl. 148, 1 (2003). 

[8] H. V. Peiris et al., Astrophys. J. Suppl. 148, 213 (2003). 

[9] D. N. Spergel et al., Astrophys. J. Suppl. 170, 377 (2007). 
[10] S. Perlmutter et al., Astrophys. J. 517, 565 (1999). 
[11] A. G. Riess et al, Astrophys. J. 116, 1009 (1998). 
[12] U. Scljak, A. Slosar, and P. McDonald, JCAP 0610, 014 (2006). 
[13] M. Tcgmark ct al., Phys. Rev. D 74, 123507 (2006). 
[14] T. Giannantonio ct al., Phys. Rev. D 74, 063520 (2006). 
[15] D. J. Eisenstein, Astrophys. J. 633, 560 (2005). 
[16] N. C. Tsamis and R. P. Woodard, Nucl. Phys. B 474, 235 (1996). 
[17] N. C. Tsamis and R. P. Woodard, Ann. Phys. (NY) 253, 1 (1997). 

[18] V. F. Mukhanov, L. R. W. Abranio, and R. H. Brandcnbcrgcr, Phys. Rev. Lett. 78, 1624 (1997). 

[19] L. R. W. Abramo, R. H. Brandcnbcrgcr, and V. F. Mukhanov, Phys. Rev. D 56, 3248 (1997). 

[20] L. R. W. Abramo and R. P. Woodard, Phys. Rev. D 60, 044010 (1999). 

[21] B. Losic and W. G. Unruh, Phys. Rev. D 72, 123510 (2005). 

[22] S. B. Giddings, D. Marolf, and J. B. Hartle, Phys. Rev. D 74, 064018 (2006). 

[23] L. R. W. Abramo and R. P. Woodard, Phys. Rev. D 65, 063515 (2002). 

[24] G. Goshnizjani and R. Brandcnbcrgcr, Phys. Rev. D 66, 123507 (2002). 

[25] J. Garriga and T. Tanaka (2007), arXiv:0706.0295 [licp-th]. 

[26] D. Espriu, T. Muhamaki, and E. C. Vagcnas, Phys. Lett. B 628, 197 (2005). 

[27] J. A. Cabrer and D. Espriu (2007), arXiv:0710.0855 [gr-qc]. 

[28] J. F. Donoghue, Phys. Rev. Lett. 72, 2996 (1994). 

[29] J. F. Donoghue, Phys. Rev. D 50, 3874 (1994). 

[30] N. E. J. Bjerrum-Bohr, J. F. Donoghue, and B. R. Holstcin, Phys. Rev. D 68, 084005 (2003), crratum-ibid. D 71, 069904 
(2005). 

[31] N. D. Birrell and P. C. W. Davies, Quantum fields in curved space (Cambridge University Press, Cambridge, 1994). 
[32] R. M. Wald, Quantum field theory in curved spacetime and black hole thermodynamics (The University of Chicago Press, 
Chicago, 1994). 

[33] E. E. Flanagan and R. M. Wald, Phys. Rev. D 54, 6233 (1996). 

[34] T. S. Bunch and P. C. W. Davies, Proc. R. Soc. London A 360, 117 (1978). 

[35] P. R. Anderson, W. Eaker, S. Habib, C. Mohna-Pan's, and E. Mottola, Phys. Rev. D 62, 124019 (2000). 

[36] L. Parker, in Recent Developments m Gravitation, edited by S. Deser and M. Levy (Plenum Press, New York, 1979). 

[37] A. Campos and E. Verdaguer, Phys. Rev. D 49, 1861 (1994). 

[38] A. Campos and E. Verdaguer, Int. J. Theor. Phys. 36, 2525 (1997). 

[39] E. Calzetta, A. Campos, and E. Verdaguer, Phys. Rev. D 56, 2163 (1997). 

[40] G. Leibbrandt, Rev. Mod. Phys. 47, 849 (1975). 

[41] R. Martin and E. Verdaguer, Phys. Lett. B 465, 113 (1999). 

[42] R. Martin and E. Verdaguer, Phys. Rev. D 60, 084008 (1999). 

[43] B. L. Hu and E. Verdaguer, Class. Quant. Grav. 20, Rl (2003). 

[44] B. L. Hu and E. Verdaguer, Living Rev. Rel. 7, 3 (2004). 

[45] E. Calzetta and B. L. Hu, Phys. Rev. D 49, 6636 (1994). 

[46] B. L. Hu and A. Matacz, Phys. Rev. D 51, 1577 (1995). 

[47] B. L. Hu and S. Sinlia, Phys. Rev. D 51, 1587 (1995). 

[48] A. Campos and E. Verdaguer, Phys. Rev. D 53, 1927 (1996). 

[49] A. Roura and E. Verdaguer, Phys. Rev. D 60, 107503 (1999). 

[50] G. T. Horowitz, Phys. Rev. D 21, 1445 (1980). 

[51] A. Campos, R. Martin, and E. Verdaguer, Phys. Rev. D 52, 4319 (1995). 
[52] L. Schwartz, Theory des distributions (Hermann, Paris, 1966). 
[53] F. Cooper and E. Mottola, Phys. Rev. D 40, 456 (1989). 
[54] J. Z. Simon, Phys. Rev. D 41, 3720 (1990). 
[55] J. Z. Simon, Phys. Rev. D 43, 3308 (1991). 



19 



J. Z. Simon, Phys. Rev. D 45, 1953 (1992). 

L. Parker and J. Z. Simon, Phys. Rev. D 47, 1339 (1993). 

B. L. Hu and A. Roura (2007), accepted for publication in Phys. Rev. D, arXiv:0708.3046 [gr-qc]. 
V. P. Frolov and A. I. Zelnikov, Phys. Rev. D 29, 1057 (1984). 

E. E. Flanagan, Class. Quant. Grav. 21, 3817 (2004). 

T. Futamase, T. Rothman, and R. Matzner, Phys. Rev. D 39, 405 (1989). 

S. Habib, C. Molina-Pan's, and E. Mottola, Phys. Rev. D 61, 024010 (2000). 

P. R. Anderson, C. Mohna-Pan's, and E. Mottola, Phys. Rev. D 67, 024026 (2003). 

P. R. Anderson, C. Molina-Paris, and E. Mottola, in preparation. 

J. F. Donoghue, in The Eighth Marcel Grossmann Meeting, edited by T. Piran and R. Rufflni (World Scientific, Singapore, 
1999), gr-qc/9712070. 

C. P. Burgess, Living Rev. Rel. 7, 5 (2004). 

S. Weinberg, Phys. Rev. D 72, 043514 (2005). 
S. Weinberg, Phys. Rev. D 74, 023508 (2006). 

B. L. Hu, A. Roura, and E. Verdaguer, Phys. Rev. D 70, 044002 (2004). 

A. Roura and E. Verdaguer, Int. J. Theor. Phys. 38, 3132 (1999). 

B. Losic and W. G. Unruh, Phys. Rev. D 74, 023511 (2006). 
J. Schwinger, J. Math. Phys. 2, 407 (1961). 

L. V. Keldysh, Zh. Eksp. Teor. Fiz 47, 1515 (1964), [Sov. Phys. JETP 20, 1018 (1965)]. 

K. Chou, Z. Su, B. Hao, and L. Yu, Phys. Rep. 118, 1 (1985). 

R. D. Jordan, Phys. Rev. D 33, 444 (1986). 

E. Calzetta and B. L. Hu, Phys. Rev. D 35, 495 (1987). 

R. P. Feynman and F. L. Vernon, Ann. Phys. (NY) 24, 118 (1963). 

R. P. Feynman and A. R. Hibbs, Quantum mechanics and path integrals (McGraw-Hill, New York, 1965). 
J. Garriga, O. Pujolas, and T. Tanaka, Nucl. Phys. B 605, 192 (2001). 
J. Garriga, O. Pujolas, and T. Tanaka, Nucl. Phys. B 655, 127 (2003). 



